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ABSTRACT 

We develop a new method which measures the projected density distribution w p (r p )n 
of photometric galaxies surrounding a set of spectroscopically-identified galaxies, and 
simultaneously the projected cross-correlation function w p (r p ) between the two popu- 
lations. In this method we are able to divide the photometric galaxies into subsam- 
ples in luminosity intervals even when redshift information is unavailable, enabling us 
to measure w p (r p )n and w p (r p ) as a function of not only the luminosity of the spec- 
troscopic galaxy, but also that of the photometric galaxy. Extensive tests show that 
our method can measure w p {r p ) in a statistically unbiased way. The accuracy of the 
measurement depends on the validity of the assumption inherent to the method that 
the foreground /background galaxies are randomly distributed and are thus uncorrelated 
with those galaxies of interest. Therefore, our method can be applied to the cases where 
foreground/background galaxies are distributed in large volumes, which is usually valid 
in real observations. 

We have applied our method to data from the Sloan Digital Sky Survey (SDSS) 
including a sample of 10 5 luminous red galaxies (LRGs) at z ~ 0.4 and a sample of 
about half a million galaxies at z ~ 0.1, both of which are cross-correlated with a deep 
photometric sample drawn from the SDSS. On large scales, the relative bias factor of 
galaxies measured from w p (r p ) at z ~ 0.4 depends on luminosity in a manner similar 
to what is found for those at z ~ 0.1, which are usually probed by autocorrelations 
of spectroscopic samples in previous studies. On scales smaller than a few Mpc and 
at both z ~ 0.4 and z ~ 0.1, the photometric galaxies of different luminosities exhibit 
similar density profiles around spectroscopic galaxies at fixed luminosity and redshift. 
This provides clear observational support for the assumption commonly-adopted in halo 
occupation distribution (HOD) models that satellite galaxies of different luminosities 
are distributed in a similar way, following the dark matter distribution within their host 
halos. 
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INTRODUCTION 



In cold dark matter dominated cosmological models, dark matter halos form in density peaks 
in the universe under the influence of gravity, and thus are clustered in a different way from the 
underlying dark matter. In other words, they are biased in spat ial distribution relative to dark 
matter (e.g. iMo fc White! Il996l: Ij ine 19981: ISeljak Warrenl 12004 ) . Galaxies are believed to form 



inside these halos (IWhite 



Rccs 



Kaiser 



1978), and thus their spatial distribution is also biased with 



1984; 



Davis et al. 



1985 



Bardeen et al.lll986l ). On large scales 



respect to dark matter (e.g. 
(> lOMpc), such biasing is nearly linear and the clustering of dark matter is well described by 
linear perturbation theory. On smaller scales, in contrast, galaxies do not trace dark matter simply. 
Complicated physical processes involved in galaxy format ion and evolution hav e to be considered 



if on e desires to fully un derstand galaxy clustering (e.g. IWhite Frenkl Il99ll : iKauffmann et al 
19991 ; IColberg et al.l l2000l ) . This leads galaxy clustering and biasing to depend on a variety of 



factors including spatial scale, redshift and galaxy properties. Therefore measuring the clustering 
of galaxies as a function of their physical properties over large ranges in spatial scale and redshift 
is helpful for understanding how galaxies have formed and evolved. 

Recent large r edshift surveys, in particular the Two Degree Field Galaxy Re dshift Survey 
(|Colless et al.l boOll . 2dFGRS) and the Sloan Digital Sky Survey ([York et al.l [joool . SDSS), have 
enabled detailed studies on galaxy clustering in the nearby universe. These studies have well 
established that the clustering of galaxies depen ds on a variety of properties, such as luminosity , 



stellar mass, color, spectral type, and morphology ( Norberg et al 



Zehavi et al.ll2002l . 12003 : iGoto et al.ll2003l : iLi et al 



2006; 



2001 



Zehavi et al 



20021 : iMadewick et al.ll2003l : 



2010j) . More luminous (mas- 



sive) galaxies are found to cluster more strongly than less luminous (massive) galaxies, with the 
luminosity (mass) dependence being more remarkable for galaxies brighter than L* (the characte r- 
istic luminosity of galaxy luminosity function described by a Schechter function. ISchechterl (|1976l )). 
Moreover, galaxies with redder colors, older stellar populations and more bulge-dominated structure 
show higher clustering amplitudes and steeper slopes in their two-point correlation functions. 

There have also been recen t studies on galax y clustering at higher redshifts. At z ~ 1, the 
DEE P2 Galaxy Redshift Survey (|Davis et alJboO.^ and the VIMOS-VLT Deep Survey \he Fevre et al 
20051 . VVDS) have shown that galaxy clustering depends on luminosity, stellar mass, c olor, spec- 



tral type and morphology, large l y consistent with what are found for t he local universe (ICoil et al 



2004 



2006 



Meneux et al 



2006 



Coil et al 



2008; 



contrast, the zCOSMOS survey (jLillv et al.H2007n shows no cle ar luminosity dependence of galaxy 



Meneux et al. 



2008 



de la Torre et al 



20091 ). In 



clustering over redshift range 0.2 ^ z ^ 1 ([Meneux et al.ll2009l ). More surprisingly, the projected 
two-point auto-correlation fu nction w p (r p ) derived from the zCOSMOS is significantly higher and 
flatter than from the VVDS dMeneux et alJbood . bood ). 



The observational measurements of galaxy clustering at both z ~ and z ~ 1 as described 



above have been widely used to test theories of galaxy formation (e.g. IKauffmann et al. 



1997 



Benson et al.1 12000 : ILi et al.1 120071 : iGuo et al.1 l2010i ) , as well as to quantify the evolution of galaxy 
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clust ering from high to low redshifts (e.g. IZheng et al.l 120071 ; iMeneux et al.l 120081 ; IWang &; Jing 



2010). Galaxy clustering has also been used to constrain halo occupation distribution(HOD) 



hence useful clues for understanding galaxy formation (e.g. 


Jing et al. 


1998; . 


fing & Boernei 


1998; 


Peacock & Smith 


2000 


Ma & Frv 


2000; 


Seliak 


2000: 


Scoccimarro et al 


2001; 


Berlind Sz Weinberg 


2002; 


Coorav & Sheth 


2002; 


Yang et al. 


2003 


; 


Zheng et al. 


2005 


: Tinker et al 


2005). 



At intermediate redshifts (0.2 < z < 1), progress on measuring galax y clustering has been 



relatively hampered by the lack of su i table data sets. A few studies (e.g. IShepherd et al.l 12001 



Carlberg et al.l l200ll ; iFirth et al.l [2002J; iPhleps et ajj 120061 ) have measured galaxy clustering as a 
function of color which are in broad agreement with results found for the local universe. However, 
the dependence of clustering on luminosity, which is well seen in the local universe, has not been fully 
established at these intermediate redshifts, very likely due to the limited size of the spectroscopic 
samples. These samples usually cover small area on the sky, suffering from both sampling noise 
and large-scale structure noise (the so-called cosmic variance effect). 

In this paper, rather than measuring the auto- correlation of these galaxies as in most pre- 
vious studies, we develop a new method for estimating the projected two-point cross- correlation 
function w p {r p ) between a given set of spectroscopically identified galaxies and a large sample of 
photometric galaxies. In brief, we first estimate the angular cross-correlation function between the 
spectroscopic and the photometric samples. We then determine the projected, average number 
density distribution w p {r p )n of the photometric galaxies surrounding the spectroscopic objects, as 
well as the projected two-point cross-correlation function w p (r p ). The photometric sample is usu- 
ally the parent sample of the spectroscopic galaxies, but goes to much fainter limiting magnitudes. 
The spectroscopic sample could be clusters (or groups) of galaxies, central galaxies of dark matter 
halos such as the luminous red galaxies (LRGs) in the SDSS, quasars, or any spectroscopic galaxy 
populations of interest. Our method can yield a measurement of the projected density distribution 
of galaxies with certain physical properties (such as luminosity, color, etc.) around spectroscopic 
objects of certain properties. In this paper we focus on presenting our methodology and limit the 
application to galaxies of different luminosities. We plan to examine the dependence of w p {r p )n 
and w p {r p ) on other properties (color, morphology, etc.) in future work. 

Previous studies of satellite galaxy distribution around relatively bright galaxies are mostly lim 



1991 



Lor rimer et al. 1994 



e.g. 


Lake & Tremaine 


1980; 


Phillipps &; Shanks 


1987 


Vader &; Sandagc 


Sales & Lambas 


2005; 


Chen et al. 


2006 


)• 


Masiedi et al. ( 


2006 


) and 



Zehavi et al.l (|2005l ) have recently investigated cross-correlations between spectroscopic and imag- 
ing galaxy samples at intermediate redshift (0.2 < z < 0.4), but with different methods and focuses 
Here we apply our method to a deep, photometric galaxy catalogue and a spectro scopic LRG sam- 
ple a t z ~ 0.4, both of which are drawn from the final data release of the SDSS (jAbazajian et al 



2009) . The LRGs are expected to be the central galaxy of their host dark matter halos. There- 
fore, by measuring w p {r p )n on scales smaller than a few Mpc, we yield an estimate of the density 
distribution of satellites galaxies within their host halo, as well as its dependence on luminosities of 
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both central and satellite galaxies. On larger scales, our analysis leads to a measurement of linear 
relative bias factor for photometric galaxies of different luminosities. 

We describe our galaxy samples in § [2] and present our methodology in § [3j Applications to 
SDSS data are presented in § [H We summarize and discuss in the last section. Throughout this 
paper we assume a cosmology with fl m = 0.3,$7a = 0.7 and Hq = 100/ikms _1 Mpc~ 1 (h = 1). 



2. Data 



2.1. The LRG sample at intermediate redshift 



The LRG sample is constructed from the SDSS data release 7 (DR7 lAbazajian et al.l 12009 ) . 
consisting of 101,658 objects with spectroscopically measured redshift in the range 0.16 < z < 0.47, 
absolute magnitude limited to —23.2 < Mo.z g < —21.2 and redshift confidence parameter greater 
tha n 0.95. Here Mo.3„ is t he g-band absolute magnitude K- and -E-corrected to redshift z = 0.3 
(see lEisenstein et al.1 l200ll for references). We further select those LRGs that are expected to be 
the central galaxy of their host dark matter ha los, using a method similar to that adopted in 
Reid fc Spergel (hood ) and lokumura et al.l (boOfll ). We use linking lengths of 0.8 /i _1 Mpc and 20 
h M.pc for separations perpendicular and parallel to the line of sight when linking galaxies into 
groups. This leads to a total of 93802 central galaxies (about 92.3% of th e initial LRG catalogue), 
covering a sky area which is almost the same as that of lKazin et al.l (|2010l ). Prom this catalogue we 



-21.8 < M . 3g < 



-21.2) 



select five samples in two luminosity intervals (—23.2 < Mo.3 g < —21.8 and 
and in three redshift intervals (0.16 < z < 0.26, 0.26 < z < 0.36 and 0.36 < z < 0.46). Details of 
our samples are listed in Table [TJ These samples so selected are volume limited, exce pt Sample L4 



which is approximately, but not perfectly volume limited as can be seen from fig. 1 of lZehavi et al 



(|2005l ). Figure Q] shows the redshift distributions of our LRGs in the two luminosity intervals. 



2.2. The low-redshift galaxy sample 

Our spectroscopic galaxy sample in the local Universe is constructed from the New York 
University Value Added Galaxy Catalog (NYU-VAGC )0, which is built by blanton et all (fcood ) 



based on the SDSS DR7. From the NYU-VAGC we select a magnitude-limited sample of 533,731 
objects with 0.001 < z < 0.5 and r-band Petrosian magnitude in the range 10.0 < r < 17.6. The 
sample has a median redshift of z = 0.09, with the majority of the galaxies at z < 0.25. The 
galaxies are divided into five non-overlapping redshift bins, ranging from z = 0.03 to z = 0.23 with 
an equal interval of Az = 0.04. The galaxies in each redshift bin are further restricted to various 
luminosity ranges, giving rise to a set of eight volume limited samples as listed in Table The r- 



4 http://sdss.physics.nyu.edu/vagc/ 
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band absolute magnitude Mo.i r is K- and i?-corrected to its value at z = 0.1 following |Blanton_et_al 



(|2003l ) (hereafter B03). Figure [2] shows the redshift distribution of the galaxies falling into the three 
luminosity ranges which are used to select the samples. These samples by construction are at lower 
redshifts when compared to the LRG samples, allowing us to make use of photometric galaxies 
for our analysis over wide ranges in luminosity and redshift. We don't attempt to select central 
galaxies as done for LRGs above, as it is not straightforward to do so. Thus one should keep in 
mind that by using the low-redshift samples selected here we will measure density profiles and 
projected correlations for general populations of galaxies, not only for central galaxies. 



2.3. The photometric galaxy sample and random samples 



We construct our photometric galaxy sample from the datasweep catalogue which is included 
as a part of the NYU-VAGC. T his is a compr e ssed v ersion of the full photometric catalogue of 
the SDSS DR7 that was used bv blanton et al.1 (|2005l ) to build the NYU-VAGC. It contains only 
decent detections and includes a subset of all photometric quantities, which is enough for our anal- 
ysis. Starting from the datasweep catalogue, we select all galaxies with r-band apparent Petrosian 
magnitudes in the range 10 < r < 21 after a correction for Galactic extinction and with point 
spread function and model fluxes satisfying f mo dei > 0.875 x /psf in ah five bands. In order to 
select unique objects in a run that are not at the edge of the field, we require the RUN PRIMARY flag 
to be set and the RUN EDGE flag not to be set. Finally we a lso require the galaxies to be located 
within target tiles of the Legacy Survey (jBlanton et al.l 12003). This procedure results in a sample 
of ~ 21.1 million galaxies. 



As shown by iRoss et al.1 (|2007l ). the datasweep catalogue needs to be properly masked, oth- 
erwise the angular correlation function obtained would be falsely flat on large scales. We describe 
the SDSS imaging geometry in terms o f disjoint spherical polygons (IHamilton &: Tegmarkl 120021 ; 
Tegmark et al.ll2002l : iBlanton et al.ll2005l ). which accompany the NYU-VAGC release and are avail- 
able from the NYU-VAGC website. We exclude all polygons (and thus the galaxies located within 
them) which contain any object with seeing greater than 1.5" or Galactic extin ction A r > .2. W e 
also exclude polygons that intersect the mask for galaxy M101 as described in lRoss et al.l (|2007l ). 
As a re sult, a fract i on of about 14.32% of the total survey area has been discarded, slightly larger 
than in lRoss et al.l (|20Q7I ) where the authors exclude image pixels rather than polygons with less 
critical criteria than adopted here. We also restrict ourselves to galaxies located in the main con- 
tiguous area of the survey in the northern Galactic cap, excluding the three survey strips in the 
southern cap (about 10 per cent of the full survey area). These restrictions result in a final sample 
of ~ 19.7 million galaxies. 

We have constructed a random sample which has exactly the same geometry and limiting 
magnitudes as the real photometric sample. This is done by generating sky positions (RA and Dec) 
at random within the polygons covering the real galaxies. In this work both the photometric sample 
and the random sample are cross-correlated with a given spectroscopic sample to estimate the 
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two-point angular cross-correlation function w(6) between the spectroscopic and the photometric 
samples. 



3. Methodology 

In our method we estimate in the first place the angular cross-correlation function w(9) between 
a given set of spectroscopic galaxies selected by luminosity and redshift (samples listed in Tables Q] 
and E]) and a set of photometric galaxies that, if they were at the redshift of the spectroscopic 
sample, would be expected to fall in a given luminosity range. Next, we convert w(9) to determine 
the projected density distribution w p (r p )n of the photometric galaxies around the spectroscopic 
galaxies, from which we further estimate the projected cross-correlation function w p {r p ) between 
the two populations. In this section we describe how we select photometric galaxies in a given 
luminosity range, followed by description of our measures of w{9) as well as the way of determining 
w p (r p )n and w p (r p ). 



3.1. Selecting photometric galaxies according to luminosity and redshift 

Considering a sample of spectroscopic galaxies with absolute magnitude and redshift in the 
ranges M Sj i < M s < M s ^ and z\ < z s < Z2, we want to measure the cross-correlation of this 
sample with a set of photometric galaxies with absolute magnitude M Pi i < M p < M p ^- Due to 
the lack of redshift information for the photometric sample, it is not straightforward to determine 
which galaxies should be selected in order to have a subset falling in the expected luminosity 
range. One can overcome this difficulty by the fact that the cross-correlation signal is dominated 
by those photometric galaxies that are at the same redshifts as the spectroscopic objects, while 
both foreground (below z\) and background (above z<i) galaxies contribute little. This is reasonably 
true when the redshift interval Z2 — z\ and the projected physical separation r p in consideration are 
substantially small, for the clustering power decreases rapidly (approximately a power law) with 
increasing separation. On large scales, projection effect due to contamination of foreground and 
background galaxies becomes relatively large (we will discuss more about this point in § 13.3. lj) . With 
this assumption in mind we restrict ourselves to photometric galaxies with apparent magnitude 
m P) i < m p < m Pi 2, the magnitude range for photometric galaxies to have absolute magnitude in 
the range M p i < M p < M Pt 2 at redshift z\ < z p < Z2, when estimating angular cross-correlation 
functions. 

When calculating the apparent magnitude for a given abs olute magnitude an d redshift, we 



have adopted the empirical formula of iT-correction presented by IWestra et al.1 (|2010l ) , which works 
at r-band as a function of observed g — r color and redshift. Since .fC-correction value changes 
slowly with redshift, we adopt (z\ + Z2)/2 as the input redshift value when applying the formula for 
simplicity. In this paper we adopt the 01 r-band luminosity function from B03, and so we convert 
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the apparent magnitude in r to that in 01 r using the analytical conversion formula provided by 
Blanton fc Roweisl fcooj ). 



3.2. Measuring angular correlation functions 

We use two estimators to measure the angular cross correlation function between a spectro- 
scopic sam ple and a photometric sample. For small separations (9 < 1000") we use the standard 
estimator (jDavis Peebleslll983l ) 

QD(9) 



w{9) 



1, 



(1) 



where 9 is the angular separation, and QD(9) and QR{9) are the cross pair counts between the 
spectroscopic sample and the photometric sample, and between the same spectroscopic sample and 
the random sample. Note that QR is normalized according to the ratio of the size of the photometric 
and rando m samples. Fo r separations larger than 9 ~ 1000", we instead use a Hamilton-like 
estimator ([Hamilton! 1 1 993l ) 

QD{9)RR{9) 



w{9) 



(2) 



QR{9)DR{9) 

where RR is the pair count of the random sample, and DR the cross pair count between the 
photometric and random samples. This estimator is expected to work better on large scales than 
the stand ard one, since it is less sensitive to uncertainties in the mean number density of photometric 
galaxies ([Hamilton! Il993l ). The two estimators differ in w{9) by 10% to 20% at 9 ~ 1000" in our 
cases. On smaller scales the two estimators give almost identical results, with difference at a few 
percent level and well within error bars. In order to reduce computation time, we apply the standard 
estimator to a random sample of 30 million points for separations 9 < 1000", while a smaller sample 
of 0.9 million random points and the Hamilton-like estimator are used for larger separations. 



3.3. Converting w{9) to w p (r p )n and w(r p ) 

Given a measurement of w(9) we estimate the corresponding projected cross-correlation func- 
tion w p (r p ) and projected density profile w p (r p )n in the following two ways. In this subsection the 
photometric and spectroscopic galaxy samples being considered are named Sample 1 and Sample 
2, respectively. 
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3.3.1. Direct conversion from w{9) to w p (r p ) 

The relation between angular correlation function w{6) and real space correlation function £(r) 
is given by (see Peebles 1980) 

ia\ Jo° xlxldxidx 2 alaln 1 n 2 ^(r 1> 2) 

w{6) = ' (3) 

where a and x stands for scale factor and comoving distance respectively; r\ 2 is the real space 
separation between Sample 1 and Sample 2 galaxies; M\ and M 2 are the surface number densities 
of the two samples; n\ and n 2 are their comoving spatial number densities. Taking n 2 as a sum of 
Dirac delta functions, we have 

n 2 = ^<5(ri -r k ), (4) 

k 

where f* k stands for galaxy positions in Sample 2. Thus Eqn. (|3|) becomes 

w{6) = mN 2 (5) 



iViiV 2 

Efc^diVx/d^K^dz! 



(6) 
(7) 



where 



n,k = yrf +rl - 2rir k cos{6), (8) 

and r\ and r k are the comoving distances for galaxies in Sample 1 and the feth galaxy in Sample 2. 
Here N\ and N 2 are the total number of objects in the two samples. Let N k denote the number of 
galaxies with approximately the same distance r k (or redshift z k ) in Sample 2, then Eqn. ([7]) can 
be written as 

_ E k N k j OO dN 1 /dz 1 gr 1 , k )dz 1 
{) ~ N\N 2 ■ [) 

If the redshift bin of Sample 2 is thin enough, all the galaxies within it can be regarded as at the 
same redshift. This gives rise to a much simplified relation between w(6) and £(r): 

w{e) = j?w/wm_ (10) 

On the other hand, the relation between w p {r p ) and £(r) is known to be 

w p (r p ) = 2 g(r) ^— = / ^r 1 , 2 )—dz 1 , (11) 

Jr v / r 2 _ r 2 J Z{ aZl 

where D is the comoving distance of Sample 1 galaxy, r p the projected physical separation, zi and 
z u the lower and upper limits of the redshift interval in consideration. If dN\/ dz\ and dD/ dz\ 



- 9 - 



change sufficiently slowly with redshift when compared to C( r ik) as a function of r% where ri& 
depends on z\, we can take dN\/ dz\ and dD/ dz\ out of the integral. Thus for a thin redshift bin 
the ratio between w(9) and w p {r p ) is simply approximated by 



w{6) _ dNx/dzx 
w p (r p ) N\ dD/ dzi 



(12) 



*1 ^med 



where r p = 2Dsin(2) and z me d is the median of the redshift range. In practice we instead calculate 
the ratio in the following way to take into account the redshift dependence (though very weak) of 
dNi/dzi and dN 2 /dz\: 



w{6) 
w p (r p ) 



L 



Zl 



djVi/dgi 



N\ AD I Az\ 



Zl = 



; dN 2 / dz) dz 



J z z ;(dN 2 /dz) dz 



(13) 



We also require the angular separation 6 to change accordingly with fixed r p , and thus in 
this way the angular separation of our measured w(9) changes with the redshift of spectroscopic 
galaxies when counting galaxy-galaxy pairs. To be more specific, the r p considered here range from 
r p ~ O.lMpc/h to r p ~ 25A'Ipc/h for LRGs, and from r p ~ O.lMpc/h to r p ~ 20Mpc/h for galaxies 
in the low-redshift sample, with 13 and 12 intervals of equal size in logarithmic space. Quantities 
in the right side of Eqn. (|13j) are determined either from data catalogue directly (N 2 and dN 2 / dz) 
or from the luminosity function analytically ( dAq / dz) . 



In order to understand to what extent Eqn. (|13p holds, we have performed two tests. In the 
first test, we c alculate a linear power spectrum Pi(k) using the CMBFAST code for the cosm ology 



adopted here (jSeljak k, Zaldarriagalll996l : IZaldarriaga et al.lll998l : Zaldarri aga k, Seliak 



2000), from 



which we calculate the nonlinear power spectrum P n i(k) following iPeacock &: Doddsl (j 19961 ). The 
real-space correlation function £(r) is then obtained by Fourier transforming P n i(k). The amplitude 
of £(r) is arbitrarily given which has no effect on the ratio of w{6) and w p (r p ). Next, the redshift 
distribution dNi/ dz± for Sample 1 galaxies in a certain magnitude range is calculated analytically 
from the luminosity function of B03. Using Eqn. ([9]) and (jlip we determine w{9) and w p {r p ) for 
this magnitude range, giving rise to the true value of their ratio which we denote as ratiotme- When 
integrating the right part of Eqn. ([9]) we fix r p and let the binning of 6 vary accordingly. We also 
calculate an approximated value for the same ratio using Eqn. (|13p . which we denote as ratio ana i y 
and compare to the true ratio in order to test the validity of Eqn. (fT3|) . 

Figure [3] shows the relative difference between the approximated and the true values of the 
w(6)/w p (r p ) ratio, (ratio ana i y — ratiotme) /ratiotme, for Sample 2 galaxies with 0.07 < z 2 < 0.078 
and —23.0 < M 2 < —21.0, and Sample 1 galaxies in several absolute magnitude intervals (as indi- 
cated in each panel). The approximated ratio agrees quite well with the true value, at 1% accuracy 
or better, for separations r p < 10 Mpc and for all luminosities considered. The discrepancy increases 
at larger separations, but well below 3% level even at the largest scale probed (~ 30Mpc). This 
discrepancy mainly comes from the distant-observer approximation adopted here. The accuracy of 
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Eqn. ()13p is expected to be better for higher redshifts where the distant-observer approximation 
works better. We thus conclude that our approximation in Eqn. (|13[) works at substantially high 
accuracies for our purpose. 

A second test that we have done is to apply our method to spectroscopic samples. Simply 
speaking, w p (r p ) for a spectroscopic sample can be measured with the redshift information. It can 
also be estimated with our method without using the redshift information for the photometric sam- 
ple. By comparing the two w p (r p ) estimates we are able to understand how well our method works. 
For simplicity we consider here a specific case in which a spectroscopic sample of given luminosity 
and redshift ranges is cross-correlated with spectroscopic (for the true w p (r p )) or "photometric" 
(for the w p {r p ) obtained by our method) galaxies in the same ranges. Thus the w p {r p ) are reduced 
to auto-correlation functions. 

The result of this test is shown in Figure HI where we plot the true w p {r p ) in blue curves and the 
approximated one in red for SDSS Main galaxies (left column) and LRGs (right column) at different 
luminosities and redshifts (indicated in each panel). For LRGs we see good agreement between the 
two measurements on all scales and at all redshifts probed (with the difference < 20%). A similar 
agreement is seen for the low-redshift samples of —20 < Mo.i r < —19 and —21 < Mo.i r < —20. All 
these results are very encouraging. 

However, there is large difference (~ 50%)between the results of the two methods on large 
scales (> 5Mpc/h) for the brightest (—22 < Mo.i r < —21) low reshift sample. The deviation may 
be caused by a coincident correlation between foreground galaxies in the photometric sample and 
the spectroscopic sample. We have performed a further analysis by estimating the cross-correlation 
with the foreground, the background, and the right redshift interval separately, for three low- 
redshift spectroscopic samples (corresponding to the left-hand panels in Fig. 2]). We find that the 
contamination comes mainly from the foreground for the brightest sample (—22 < Mo,i r < —21), 
and from the background for the faintest sample (—20 < Mo.i r < —19). For both samples, the 
projected cross-correlation w p (r p ) with the foreground (for the brightest sample) or the background 
(for the faintest sample) shows weak dependence on scale. When compared to the true w p (r p ), the 
cross-correlation with the foreground/background is negligible on small scales, ~ 50% smaller at 
~ 10Mpc//i and compatible at ~ 20Mpc//i. This result clearly shows that the clustering pattern 
of the forground /background can contaminate the angular cross correlation function stochastically. 
This also explains why we can measure the projected correlation function for the LRG sample 
accurately, because the foreground/background galaxies are in big cosmic volumes and thus have 
weak correlations themselves. 

We conclude that the accuracy of our method relies on the key assumption that foreground/background 
galaxies have weak correlation with the spectroscopic galaxies. This assumption is valid for many 
real observations, especially for spectroscopic samples at intermediate or high redshift. This is why 
we can recover the project correlation function for LRG samples on all scales. For the low red- 
shift samples, our method works for the sample of luminosity M*, since the forground/background 
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galaxies are relatively small in number compared with those at the redshift of the spectroscopic 
sample. For the bright low redshift sample, the foreground galaxies are located in a small volume 
and are more numerous, and their clustering pantern can bias the estimation of the projected func- 
tion. This effect is smaller for small scales, which is the reason why we can measure the projected 
function accurately on scales smaller than ~ 1 Mpc//i. 



3.3.2. Indirect conversion through w p {r p )n 

We propose a second method here for estimating w p {r p ). Rather than directly converting w{6) 
to w p (r p ), we first convert the former to a projected density profile w p (r p )ni, from which we then 
estimate w p {r p ) by calculating analytically the spatial number density n\ from the luminosity func- 
tion. w p (r p )n\ is obtained as a byproduct without suffering from uncertainties in galaxy luminosity 
function. 

In this method we estimate a weighted angular correlation function w(9) we i g ht instead of the 
traditional function w(9) as discussed above. This is measured using the same estimators given in 
Eqn. (prj and (J2]), except that each spectroscopic galaxy in Sample 2 is weighted by D~ 2 , inverse 
of the square of comoving distance for Sample 2 galaxies. It can be easily proved with Eqn. ( [T3|) 
that for a thin redshift interval of the spectroscopic sample (Sample 2), Niw(9) we i g ht/Q equals to 
the projected density profile w p (r p )ni, i.e., 



where N\ and n\ are the number and number density of photometric galaxies in Sample 1, and 
O the total sky coverage of Sample 1. Given the projected density profile w p {r p )ni estimated 
by Eqn. (|14j) as well as the spatial number density of galaxies n\ analytically calculated from 
the luminosity function, we finally estimate the projected correlation function w p (r p ) by dividing 
w p (r p )ni by n\. We emphasize here the projected density profile w p (r p )ni does not suffer from 
uncertainties in luminosity function, because all quantities in the right side of Eqn. (|14p can be 
obtained from data. 

To calculate n\, we adopt the luminosity evolution model and the luminosity function at z = 0.1 
from B03 when doing calculation for spectroscopic galaxies at low redshifts (samples selected from 
the SDSS Main galaxy catalogue as listed in Table [2]). Considering that the evolution model of B03 
is based on low-redshift data (z < . 25), w hich might not be suitable for higher redshifts, we adopt 



the evolution model of iFaber et al.l (|2007l ) (here after F07) for our LRG samples. Moreover, we 
need to convert the F07 model from 5-band to the 01 r-band at which our galaxies are observed. 
Assuming that the slope of the luminosity evolution doesn't depend on waveband, we obtain the 
0,1 r-band luminosity evolution model by simply shifting the amplitude of the -B-band model from 
F07 so as to have an amplitude at z = 0.1 which is equal to the amplitude of the °' 1 r-band model 
of B03. In this manner the slopes of the two models remain unmodified, which are Q = —1.23 for 
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F07 and -1.62 for B03 respectively (see their papers for details 

In conclusion, the projected cross-correlation function w p (r p ) can be measured either from 
Eqn. (|13p by direct conversion of w(9), or from Eqn. (|14p by indirect conversion through esti- 
mating w p {r p )n. After having made extensive comparisons, we found that the two methods give 
rise to almost identical results. In what follows we choose to use the second method only, as it 
simultaneously provides both w p {r p ) and w p (r p )n. 



3.4. Division and combination of redshift subsamples 

In this subsection we address an important issue which we have ignored so far. As mentioned 
above, our method for selecting photometric galaxies according to luminosity and redshift is valid 
only when the redshift interval Z2 — Z\ is small enough. However, the redshift intervals used to 
select our spectroscopic galaxy samples apparently do not satisfy this condition. For example, for 
a redshift range of 0.16 < z < 0.26 and an absolute magnitude interval of M p i — M Pt i = 0.5, the 
photometric galaxies selected will cover a much broader apparent magnitude range, m P) 2 — m p,i 
1 .7. 

Our solution is to further divide the galaxies in a given spectroscopic sample into a number of 
subsamples which are equally spaced in redshift (hereafter called redshift sub-shells) . See Tables Q] 
and [2] for the number of redshift sub-shells adopted for our samples. For a given sample, we 
measure the weighted angular cross-correlation function w(6) we i g ht (see above) for each sub-shell 
separately by cross-correlating with galaxies selected from the photometric catalogue in the way 
described above according to the expected luminosity range and the redshift range of the sub- 
shell. Each w{9) we i g ht measurement is then converted to give the corresponding projected density 
profile w p (r p )n as well as the projected cross-correlation function w p (r p ), using the second method 
described above. Estimates of these quantities for the sub-shells are then averaged to give the 
estimates for their parent sample as a whole. In this procedure each sub-shell is weighted by Vi/af, 
with VI being the comoving volume covered by the ith sub-shell and of the variance of w p {r p )n or 
w p (r p ) of the sub-shell. In order to estimate <7j we have generated 100 bootstrap samples for each 
sub-shell. The variance <7j of a sub-shell is then estimated by the la scatter between all its bootstrap 
samples. This weighting scheme ensures the averaged w p (r p )n or w p (r p ) to be determined largely 
by sub-shells with relatively large volume and high signal-to-noise ratio (S/N) measurements, thus 
effectively reducing the overall sampling noise and cosmic variance. 

In order to increase the accuracy of our method, one may want to increase the number of sub- 
shells for a given redshift range, at the cost of increasing both the sampling noise and the large-scale 
structure noise (the cosmic variance). In practice, we split a spectroscopic sample into redshift sub- 
shells by requiring Az/z < 0.1, where Az is the thickness of the sub-shells and z is the mean redshift 



5 We have repeated our analysis for LRGs, adopting Q — 



— 1.62 instead of Q = 



— 1.23, and obtained similar results. 
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of the sample. With this restriction the difference between m<i — m\ and M2 — M\ ranges from 
~ 0.25 occurring for low-redshift sub-shells to ~ 0.1 for high-redshift ones. Extensive tests show 
that our results are robust to reasonable change of the thickness of sub-shells. For instance, taking 
Sample LI from Table Q] as the spectroscopic sample, the cross-correlation function measured by 
10 sub-shells differs from the one of 5 sub-shells by at most 10% for photometric galaxies with 
-22 < M < -21.5, and by only about 3% for the faintest luminosity bin (-19.5 < M < -19.0). 



3.5. Error estimation 

We estimate the error in the averaged w p (r p )n or w p (r p ) by 



° 2 *& if ^i:> 1 

,.2 



(15) 



with 

- 2 = >;-r / 1 >;-7i * no) 




X /dof. = — ) Axi- x avg ) cr- , (17) 



where the sum goes over all the sub-shells of a given spectroscopic sample; Xi is the measurement 
of w p (r p )n or w p (r p ) of the ith sub-shell and x avg the average measurement for all the sub-shells 
as a whole; N su b is the number of sub-shells. Overall, Eqn. (fT5l) should be able to include both the 
volume effect (through factor Vi) and the sampling noise (through at), thus providing a reasonable 



estimate of the errors in our measurements. By weighting the error by y ^j-, we mean to take into 
account the large variation from sub-shell to sub-shell in some cases. 

To better understand the error contribution from different redshift sub-shells, we plot in Fig- 
ure [5] the w p (r p ) measurements for Sample LI listed in Table [TJ Different panels correspond to 
photometric galaxies in different luminosity intervals. In each panel, we plot w p {r p ) for all the five 
sub-shells with their redshift ranges indicated in the bottom-right panel. Error bars on the w p (r p ) 
curves are estimated using the bootstrap resampling technique, i.e. Oi in Eqn. (fT6|) and (fT7|) . We 
see that, for photometric galaxies at fixed luminosity, w p (r p ) measurements of different sub-shells 
are almost on top of each other, indicating that the scatter between sub-shells is fairly small. This 
again shows that the correlation functions measured with our method are insensitive to the number 
of redshift sub-shells. 

We note that the overall error increases rapidly with luminosity at the bright end. This 
reflects not only the sampling noise of the small samples, but more importantly, also an effect of 
a huge foreground population which significantly suppresses the angular cross-correlation signals. 
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Letting w'(8) be the angular correlation function between a spectroscopic sample with Z\ < z < Z2 
and a photometric sample including galaxies of all redshifts, and w{6) the one between the same 
spectroscopic sample and a sample of photometric galaxies within z\ < z < Z2, one can easily show 
that 

W ' {6) = J^ w{e) > (18) 
where Nq and Nqs are respectively the number of photometric galaxies in the full sample and in 
the redshift range z\ < z < Z2- In this case the estimated correlation signal is suppressed by a 
factor of -jyjr, a large effect in particular when Nq ^> Nqs as in bright samples. 

This is explained more clearly in Figure [6] where we plot the redshift distribution as calculated 
using the luminosity function of B03 for photometric galaxies which are selected according to 
luminosity and redshift using the method described in § 13.11 Plotted in different lines are the 
distributions for different luminosity ranges with redshift range fixed to 0.2 < z < 0.22. According 
to our method of dividing photometric sample into luminosity subsamples, photometric galaxies 
selected to serve our purpose of a certain luminosity bin have the desired luminosity at the chosen 
redshift. As can be seen from the figure, in the brightest luminosity interval (—22.5 < Mo.i r < 
—22.0, solid line) only a small fraction of the galaxies are located within the expected redshift 
range, where the majority of galaxies are intrinsically fainter. It is this population that suppresses 
the angular cross-correlations that we measure at the bright end, leading to large uncertainties in 
w p (r p ), as seen in Figure [5j 

Before we apply our method to SDSS data in the next section, we should point out that, 
although we have carefully considered both the sampling noise and the large-scale structure noise, 
our error estimation doesn't include the projection effect caused by background and foreground 
galaxies. Thus when interpreting our w p (r p )n and w p (r p ) presented below, one should keep in mind 
that their errors are underestimated to varying degrees, depending on the redshift and luminosity 
we consider. 



4. Applications to SDSS galaxies 

4.1. Projected cross-correlations and density profiles of LRGs 

In Figure [7] we show the projected density profile w p (r p )n for LRGs in different intervals of 
luminosity and redshift, as traced by surrounding photometric galaxies of different luminosities. 
Results are plotted for LRGs with —23.2 < Mo.s g < —21.8 in panels on the left and for those with 
—21.8 < Mo.ig < —21.2 on the right, with panels from top to bottom corresponding to different 
redshift bins. The w p (r p )n traced by photometric galaxies in different luminosity ranges are shown 
using different lines as indicated in the bottom-right panel. As can be seen, the projected number 
density of galaxies around central LRGs decreases as their luminosity increases. This is true for 
all redshifts and all scales probed. In particular, such luminosity dependence is weak for galaxies 
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fainter than the characteristic luminosity of the Schechter luminosity function (Mo.i r = —20.44 for 
the SDSS), and becomes remarkable for brighter galaxies. The change of the amplitude mainly 
reflects the number density of galaxies as a function of their luminosity. From the figure one can 
easily read out the number of galaxies in the host halo of the LRGs. For example, on average there 
are about 3 galaxies of —20.5 < M < —20 in the host halo of a LRG in Sample LI, assuming the 
host halo radius is 0.3/i _1 Mpc. 

In Figure [8] we show projected cross-correlation function w p {r p ) obtained from the w p (r p )n 
measurements shown in Figure for the same set of LRG samples and the same intervals of 
photometric galaxy luminosity. At fixed scale and redshift, the amplitude of w p (r p ) increases with 
increasing luminosity, a trend which has already been well established by previous studies. It is 
interesting to see from both figures that, although both w p (r p )n and w p (r p ) show systematic trends 
with luminosity in amplitude, their slope remains fairly universal for given redshift and central 
galaxy luminosity, regardless of the luminosity of surrounding galaxies we consider. This provides 
direct observational evidence that satellite galaxies of different luminosities follow the distribution 
of dark matter in the same way within their dark matter halos, an assumption adopted in many 
previous studies on HOD modeling of galaxy distribution. 

In Figured] we plot the w p (r p ) again in order to explore the evolution with redshift. Measure- 
ments of different photometric galaxy luminosities are shown in different panels, while in each panel 
we compare w p (r p ) measured at different redshifts for fixed luminosity. Note that in this figure 
we have considered the luminosity evolution of galaxies, as we aim to study how the projected 
density distribution and cross-correlation around the LRGs have evolved over the redshift range 
probed. We do not include S-correction in Figures 15171 and El because we want the results there 
to be less affected by possible uncertainty in the luminosity evolution model. From Figure we 
see significant increase in the amplitude of w p (r p ) as redshift goes from z = 0.4 to 0.2. Taking the 
—21.0 < M < —20.5 bin in the left column for example, on average the w p (r p ) amplitude differs 
by a factor of about 2 between the result of z ~ 0.2 (blue curve) and z ~ 0.4 (green curve) at 
separations r p < 0.3fo _1 Mpc, the typical boundary of LRG host halos. 

It is important to understand whether the significant evolution seen above could be explained 
purely by evolution in dark matter distribution, or additional processes related to galaxies them- 
selves are necessary. To the end we have done a simple calculation as follows. We assume that 
there is no merger occurring between galaxies or between their host halos. In this case all galaxies 
and halos have been evolving in a passive manner, and the total number of each keeps unchanged 
during the period in consideration. Thus galaxies at different redshifts are the same population 
which are hosted by the same set of halos. We calculate a dark matter density profile averaged 
over all dark matter halos that are expected to host LRGs of given redshift and luminosity, by 

Jtt:: P(r,M)n{M)&M 

Pavg{r) = ~Tm ; " , (19) 



where n(M) is the halo mass function from ISheth fc Tormenl (|1999i ). and p{r,M) is the density 
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profile of halos of mass M assumed to be in the NFW form (INavarro et al.lll997l ). We determine 
the concentration parameter c of halos following IZhao et al.l (|2009i ) . The lower and upper limits of 
halo mass (M m j n and M max ) are determined by matching the abundance of LRGs in our sampl e 
with that of dark matter halos given by the Sheth-Tormen mass function (jSheth fc Tormenl ll999), 
For this we have adopted the plausible assumption that the luminosity of a central galaxy is an 
increasing function of the mass of its halo. 

We consider three redshifts which are z =0.21, 0.31 and 0.41, approximately the median 
redshifts of our LRG samples. The dark matter density within O.S/i^ 1 Mpc around LRGs of 
-23.2 < Mo.3 g < -21.8 is predicted to increase by 24.6% from z ~ 0.4 to ~ 0.2, and 12.1% from 
z ~ 0.3 and ~ 0.2. The factor is 10.6% for -21.8 < M .3 g < -21.2 from z ~ 0.3 to ~ 0.2. It is 
clear that these predictions are much smaller when compared to what we have obtained from the 
SDSS data. As can be seen from Figure at r p < 0.3/i _1 Mpc, the clustering amplitude changes 
by a factor of ~2 between z ~ 0.4 and ~ 0.2, and by 20 - 50% between z ~ 0.3 and ~ 0.2. 

This simple calculation seems to suggest that the evolution of galaxy clustering observed in this 
work is caused not only by the evolution of the underlying dark matter, but also by the evolution 
of galaxies themselves. However, this argument should not be overemphasized. As selected in r- 
band, our photometric sample may be biased to bluer galaxies as one goes to higher redshifts. This 
selection effect must be properly taken into account when one addresses the evolution of galaxy 
clustering, but such a detailed modeling is out of the scope of our current paper. 

Figure [10] shows the relative bias factor of photometric galaxies with respect to L* galaxies 
as a function of luminosity. Taking a LRG sample, the bias factor for a sample of photometric 
galaxies at given luminosity and redshift is calculated from the amplitude of the w p (r p ) between 
the LRG and the photometric samples, normalized by the w p (r p ) of the same LRG sample with 
a photometric sample selected by —21 < Mo.i r < —20 (K— and E— corrected to z = 0.1) [§, and 
averaged over separations between r p = 2.5/i~ 1 Mpc and 10/i _1 Mpc. The bias factor so obtained 
should be virtually identical to the one estimated from the auto-correlation function of the same 
set of photometric galaxies. In the figure, curves in different colors refer to results from different 
LRG samples in Table HJ Since the photometric sample becomes somewhat incomplete for L < L* 
in the redshift range 0.36 < z < 0.46, the bias factor for this sample (the green line in the figure) 
is normalized with respect to a —21.0 < Mo.i r < —20.5 samp le instead of the —21 < Mo.i r < —20 
one. Plotted in black triangles is t he result oflLi et al.1 ( 20061 ) . Black dashed line is a fit obtained 



from the SDSS power spectrum bv lTegmark et al.l (J2004J). Relative bias factors at 0.16 < z < 0.26 
(Samples LI and L2) are well consistent with those previous studies at all luminosities, exc ept the 
bright end where our measurement is slightly lower than that from iTegmark et al.l (120041 ) . Our 
bias factor measurements show that the luminosity dependence of galaxy clustering observed in the 
local Universe is very similar to that at intermediate redshift z ~ 0.4. 



6 The absolute magnitude range for selecting L* samples varies from sub-shell to sub-shell for the photometric 
galaxies, so that the corresponding absolute magnitude range at z = 0.1 is always —21 < Mo.i r < —20. 
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4.2. Projected density profiles and clustering of low-z galaxies 



We have also measured w p (r p )n and w p (r p ) for the eight volume limited samples of low-redshift 
galaxies listed in Table [2] and for different intervals of photometric galaxy luminosity. Since the 
volume covered by these spectroscopic samples is small due to their small redshift intervals, the 



w 



PV V 



In and w p {r p ) measurements are more noisy than presented above for the LRG samples. In 
order to improve the S/N of our measurements, for a given luminosity interval of photometric 
galaxies we combine the measurements for spectroscopic samples that share a same luminosity 
interval but span different redshift ranges. When doing the combination we weight each sample by 
its comoving volume divided by the variance of the measurement, in the same way as above when 
combining redshift sub-shells. The combined w p (r p )n are plotted in Figure [TT1 We do not include 
independent plots for w p (r p ), which show behaviors very similar to w p (r p )n. 

In Figure [121 we show the corresponding relative bias factors, based on data points over 1.9 < 
r p < 10/i _1 Mpc. Results are plotted in blue, red and green curves for the three luminosity intervals 
of spectroscopic gala xies. For compariso n we also repeat the bias factors from previous work by 
Tegmark et al.1 ([2004) and lLi et al.l (12006) which are based on spectroscopic samples similar to those 
used here. Our measurements are roughly in agreement with these previous determinations. Our 
bias factors from the three samples agree with each other at the intermediate luminosities around 
L*, while showing obvious deviations at the bright and faint ends. Again, these differences should 
not be regarded as significant due to large uncertainties. 

Similar to what is found for LRGs at z ~ 0.4, the density profiles around galaxies in the local 
universe also shows quite similar slope, independent of the luminosity of surrounding photometric 
galaxies. Unlike in the LRG samples, the spectroscopic galaxies in our low-redshift samples could 
be either centrals or satellites of their host halos. However, the fraction of satellites should be 



small as the spectroscopic objects in the SDSS are relatively bright (jZheng et al.l 120071 ). Thus our 
conclusion made above for z ~ 0.4 more or less holds for the local Universe, that is, satellite galaxies 
of different luminosities are distributed within their halos in a similar way, i f the halos host central 
galaxies of similar lumi nosity. This is consistent with previous studies (e.g. IVale Ostrikeii 12004! . 



2006 



Zheng et al.ll2007l ) which have revealed a tight relation between central galaxy luminosity (L c ) 



and halo mass. Moreover, halo occupation distribution models usually assume galaxy distribution 
inside halos to trace their dark matter, based on studies of satellite distributions in simulations (e.g. 



Nagai &; Kravtsov 



2005 



Maccib et al.1 120061 ). Again, our results provide additional, clear evidence 



for supporting this assumption. 



5. Discussion and Summary 

Previous studies on galaxy clustering as a function of luminosity usually make use of spec- 
troscopic galaxy catalogue, thus are limited to relatively bright galaxies and low redshifts. In 
this work we have developed a new method which measures simultaneously the projected number 
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density profile w p {r p )n and the projected cross-correlation function w p (r p ) of a set of photomet- 
ric galaxies, surrounding a set of spectroscopic galaxies. We are able to divide the photometric 
galaxies by luminosity even when redshift information is unavailable. This enables us to measure 
w p (r p )n and w p (r p ) as a function of not only the luminosity of the spectroscopic galaxies, but also 
that of the surrounding photometric galaxies. Since photometric samples are usually much larger 
and fainter than spectroscopic ones, with our method one can explore the clustering of galaxies to 
fainter luminosities at high redshift. 

We have applied our method to the SDSS data including a sample of 10 5 luminous red galaxy 
(LRGs) at z ~ 0.4 and a sample of about half a million galaxies at z ~ 0.1. Both are cross-correlated 
with an SDSS photometric sample consisting of about 20 million galaxies down to r = 21. We have 
investigated the dependence of w p {r p )n and w p (r p ) on galaxy luminosity and redshift, by dividing 
both spectroscopic and photometric galaxies into various luminosity intervals and different redshift 
ranges. 

The conclusions of this paper can be summarized as follows. 

• We develop a new method which measures the projected density distribution w p (r p )n of 
photometric galaxies surrounding a set of spectroscopically-identified galaxies, and simulta- 
neously the projected cross-correlation function w p (r p ) between the two populations. In this 
method we are able to divide the photometric galaxies into subsamples in luminosity intervals 
even when redshift information is unavailable, enabling us to measure w p (r p )n and w p (r p ) as 
a function of not only the luminosity of the spectroscopic galaxy, but also that of the pho- 
tometric galaxy. Extensive tests show that our method can measure w p (r p ) in a statistically 
unbiased way. The accuracy of the measurement depends on the validity of the assumption 
inherent to the method that the foreground/background galaxies are randomly distributed 
and are thus uncorrelated with those galaxies of interest. Therefore, our method can be 
applied to the cases where foreground/baground galaxies are distributed in large volumes, 
which is usually valid in real observations. 

• We find that, for a spectroscopic sample at given luminosity and redshift, the projected cross- 
correlation function and projected density profile as traced by photometric galaxies show quite 
similar slope to each other, independent of the luminosity of the photometric galaxies. This 
indicates that satellite galaxies of different luminosities are distributed in a similar way within 
their host dark matter halos. This is true not only for LRGs at intermediate redshifts which 
are mostly central galaxies of their halos, but also for the general population of galaxies in 
the local Universe. Our result provides observational support for the assumption commonly- 
adopted in halo occupation distribution models that the distribution of galaxies follows the 
dark matter distribution within their halos. 

• The relative bias factors are estimated for photometric galaxies as a function of luminosity 
and redshift. In particular, we measured the bias factors of such kind, for the first time, for 
galaxies at intermediate redshift (z ~ 0.4) over a wide range in luminosity (0.3L* < L < 5L*). 
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There have been previous studies of measuring galaxy clustering by cross-correlations with 
imaging data. Although the methods and purposes of these studies are different from ours, i t 
is worthy of mentioning them and pointing out the findings in common. Eisenstein et al.l ((2005) 
measured the mean overdensity around LRGs in an earlier SDSS release by cross-correla ting 32,000 
LRG s with 16 million photometric galaxies, using the deprojecting method developed in (jEisenstein 
2003). The authors were aimed at understanding the scale- and luminosity-dependence of the 
clustering of LRGs, thus using only L* galaxies from the photometric sample as a tracer of the 
surrounding distribution. In our work we consider luminosity dependence for both LRGs and 
photometric galaxies, and this is why we have made considerable efforts on selecting photometric 
galaxies of specific luminosities. Our measurements from cross-correlation with L* galaxies show 
strong dependence in w p (r p ) amplitude on LRG luminosity, as well as an obvious transition at 
around 1 Mpc//i, which are in agreement with what those authors find. 



Masjedi et all (|2006l ) measured the cross-correlation between ~25,000 LRGs and an imaging 



sample, in order to correct for the effect of fiber collisions on their small-scale measurements of 
LRG auto-correlations. They find that the real-space auto-correlation function of LRGs, is 
surprisingly close to a r -2 power law over more than 4 orders of magnitude in separation r, down to 
r ~ 15 kpc/h. We don't have a measurement of auto-correlations £(r) or w p {r p ) for LRGs, but our 
results are not inconsistent with theirs in the sense that the projected cross-correlation of our LRG 
samples contin uously increa s es at such small scales, with a slope similar to that on large scales. In 
a recent work, I White et al.l (|201ll ) studied the clustering of massive galaxies at z ~ 0.5 using the 
first semester of BOSS data. The authors computed the projected cross-correlation between their 
imaging catalogue and the spectroscopic one, as an additional analysis to emphasize that there is 
a significant power on scales below 0.3 Mpc//i where their spectroscopy-based measurements suffer 
from fiber collisions. This is obviously consistent with what we have seen from our measurements. 

Our method can be applied and extended to many important statistical studies of galaxy for- 
mation and evolution. By dividing galaxies into red and blue populations in the photometric sample 
according to their colors, we can quantify the evolution of the blue fraction of galaxies in clusters 



and g roups from redshift 0. 4 to the present day (i.e the Butcher- Oemler effect, e. g. iGoto 



20031 : iDe Pronris et al 



et a 



2004 ). Wide deep photometry surveys, such as Pan-Starrs ( Kaiserl 2004 ) 



and LSST (jTysonl 120020 . will be available in the next years. Combining such surveys with large 
spectroscopic samples, such as BOSS LRG samples, will allow one to explo re the clustering o f 
galaxies from z = up to z = 1 for a wide range of luminosities. The WISE ([Wright et al.l 12010) 
will produce an all-sky catalogue of infrared galaxies. By combining this survey with SDSS spec- 
troscopic samples, one can study how infrared galaxies are distributed relative to the network of 
optical galaxies. 
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Table 1. LRG samples selected according to luminosity and redshift 



Sub-shell information 



Sample 


M0.3g 


z 


Num. of galaxies 


Num. of sub-shells 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


LI 


(-23.2,-21.8) 


(0.16,0.26) 


5592 


5 


0.02 


L2 


(-21.8,-21.2) 


(0.16,0.26) 


16832 


5 


0.02 


L3 


(-23.2,-21.8) 


(0.26,0.36) 


11647 


5 


0.02 


L4 


(-21.8,-21.2) 


(0.26,0.36) 


26301 


5 


0.02 


L5 


(-23.2,-21.8) 


(0.36,0.46) 


17502 


5 


0.02 



Table 2. Low-redshift volume- limited samples of galaxies selected by luminosity and redshift 



Sub-shell information 



Sample 


Mo.i r 


z 


Num. of galaxies 


Num. of sub-shells 


&z sub 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


VL1 


(-21.0,-20.2) 


(0.03,0.07) 


20553 


10 


0.004 


VL2 


(-23.0,-21.0) 


(0.03,0.07) 


5879 


10 


0.004 


VL3 


(-21.0 -20.2) 


(0.07,0.11) 


63694 


5 


0.008 


VL4 


(-23.0 -21.0) 


(0.07,0.11) 


19138 


5 


0.008 


VL5 


(-23.0,-21.0) 


(0.11,0.15) 


40226 


4 


0.01 


VL6 


(-23.0,-21.9) 


(0.11,0.15) 


2330 


4 


0.01 


VL7 


(-23.0,-21.9) 


(0.15,0.19) 


3935 


4 


0.01 


VL8 


(-23.0,-21.9) 


(0.19,0.23) 


5268 


4 


0.01 
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Fig. 1. — Redshift histograms for luminous red galaxies (LRGs) in our sample falling in the two 
luminosity intervals, as indicated, which we use to select our subsamples in Table 1. The g-band 
absolute magnitude Mo.3 g is K— and E— corrected to its value at z = 0.3. 
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Fig. 2. — Redshift histograms for spectroscopic galaxies in our low-redshift sample in three different 
luminosity ranges as indicated, which we use to select our subsamples in Table 2. The r-band 
absolute magnitude Mo.i g is K— and E— corrected to its value at z = 0.1. 
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Fig. 3. — Relative difference in w(9)/w p (r p ) ratio between the value from Eqn. (|13p and the value 
from theoretical calculation (see § 3.3.1 for details), for spectroscopic galaxies at 0.07 < z < 0.078 
and photometric galaxies at different luminosities, as indicated in each panel. 
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Fig. 4. — Plotted in each panel (red symbols connected by a red line) is w p (r p ) of spectroscopic 
galaxies at given luminosity and redshift ranges (as indicated), estimated from cross-correlation 
with photometric galaxies in the same luminosity and redshift ranges. This is compared to the 
projected auto-correlation correlation function of the same set of spectroscopic galaxies, as plotted 
in blue symbols/lines. Panels in the left-hand columns are for galaxies in the low-redshift sample 
and those in the right are for LRGs. 
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Fig. 5. — Projected cross-correlation function w p (r p ) for redshift sub-shells of Sample LI in Table[TJ 
estimated by cross correlating each sub-shell with photometric galaxies in different luminosity 
intervals as indicated in each panel. Different lines are for different sub-shells with their redshift 
ranges indicated in the bottom-right panel. 
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Fig. 6. — Redshift distribution as predicted by the luminosity function from iBlanton et al.1 ((2003) 
for photometric galaxies that are expected to fall in the indicated luminosity intervals if they were 
located in the redshift range 0.2 < z < 0.22. The two vertical lines mark the redshift range. 
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Fig. 7. — Projected density profile w p (r p )n in units of Mpc~ 2 h 2 surrounding LRGs with different 
luminosities (indicated above the figure) and redshifts (indicated in each panel), as traced by 
galaxies of different luminosities (shown in different lines in each panel and indicated in the bottom- 
right panel). 
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Fig. 8. — Projected cross-correlation function w p {r p ) measured for the same set of LRG samples 
and the same intervals of photometric galaxy luminosity as in the previous figure. 
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Fig. 9. — Each panel compares w p (r p ) measured at different redshifts (as indicated), but for fixed 
luminosity ranges for photometric and spectroscopic galaxies. The luminosity ranges of spectro- 
scopic galaxies are indicated above the figure, while those of photometric galaxies are indicated in 
each panel. 
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Fig. 10. — Relative bias factor with respect to L* as a function of luminosity. Bias factors are calcu- 
lated from the amplitude of w p (r p ) averaged over 2.5/i _1 Mpc < r < 10.0/i _1 Afpc, and normalized 
by w p (r p ) between the same LRG sample and the photometric sample with —21 < Mo.i r < —20. 
Note that for the 0.36 < z < 0.46 bin, the photometric catalogue becomes incomplete when 
M > —20.5, and so we calculate the relative bias factor with respect to the photometric sample of 
—21.0 < Mo.i r < —20.5. Different curves are for results obtained with different LRG sub samples, 



as in dicated. Black triangles and the black dashed line show previous determinations from lLi et al 



(120061 ) and (|Tegmark et al.ll2004l ). which are based on auto-correlation function or power spectrum 
of the SDSS Main galaxy sample with a mean redshift of z ~ 0.1. 
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Fig. 11. — Projected density profile w p (r p )n between spectroscopic galaxies of different luminosities 
(indicated above each column) and photometric galaxies of different luminosities and redshifts (both 
indicated in each panel). A solid black line is repeated on every panel to guide the eye. 
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Fig. 11. — Continued.. 
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Fig. 12. — Relative bias factor as a function of luminosity, measured for low-redshift samples. 
Symbols and lines are similar as in Figure 10. 



